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^ Abstract. The decay of turbulent and laminar oblique bands in the lower transitional range 

^^ of plane Couette flow is studied by means of direct numerical simulations of the Navier-Stokes 

r-^ equations. We consider systems that are extended enough for several bands to exist, thanks 

fS| to mild wall-normal under-resolution considered as a consistent and well-validated modelling 

strategy. We point out a two-stage process involving the rupture of a band followed by a 
slow regression of the fragments left. Previous approaches to turbulence decay in wall-bounded 
flows making use of the chaotic transient paradigm are reinterpreted within a spatiotemporal 



w^ perspective in terms of large deviations of an underlying stochastic process. 

i 

O 1, Context 

^^ Difficulties in understanding the transition to turbulence in wall-bounded flow arise from its 

^^ hysteretic character linked to the fact that the nontrivial turbulent regime stands at a distance 

from the laminar regime and coexists with it in some range of Reynolds numbers Re. In plane 
Couette flow - the flow between counter-translating parallel plane plates (speed ±U, gap 2/i, 
PCF for short) - the laminar (trivial) base state is linearly stable for all Reynolds numbers 
^ Re = Uh/jy and the turbulent state is uniform above Ret ^ 410. Below Rcg ^ 325, laminar flow 

"^ prevails in the long-time limit. Between Rcg and Ret, the nontrivial state displays a pattern of 

alternately turbulent and laminar oblique bands, at least in systems with in-plane dimensions 



very large when compared to the gap (Prigent et a/., 2003). The decay of turbulence in the 



vicinity of Rcg has mostly been theoretically analysed within the framework of low dimensional 
dynamical systems, using periodic domains called minimal flow units (MFUs) with periodic 



boundary conditions at distances of the order of the gap between the plates ( | Waleffe , 1997). 



This approach helps us to understand the general mechanisms sustaining the turbulent regime 
and permits us to use the concepts and tools of deterministic chaos theory. Within this 
framework, decay is described in terms of chaotic transients characterised by lifetimes with 



C^ exponentially decreasing probability distributions (Eckhardt et a/., 2008). However, from a 

practical point of view, the extended geometries of interest in experiments rather appeal to 
concepts of pattern formation and spatiotemporal chaos. Here we focus on the lower end of the 
transitional range around Rcg , on how the bands become fragmented into turbulent patches 



that finally decay (iManneville , 2011), and on the extent to which this can change our views on 



turbulence breakdown, in contrast with what was put forward at the MFU scale. 
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Figure 1. Variation of the time- average of the distance A to laminar flow as a function of R 



for A ^^ = 21 with (L^^.L^) = (110,84) and (N^^.N^) = (440,336), see (Mannevihe fc Rolland, 



2011). Sustained bands are conspicuous between Rcg below which they inevitably decay in 



the long term, and Ret above which uniform turbulence is obtained. Similar experiments for 
different in-plane resolutions and domain sizes helped us to build the diagram on the right 
displaying thresholds Rg and Rt as functions of the number Ny of Chebyshev polynomials. 
Marked by stars, the values obtained for N y = 33 (Duguet et a/., 2010) cannot be distinguished 



from experimental flndings Prigent et al. (2003). Below A^^ = 11, no bands are found like in 



Manneville (2009) resting on the low order Galerkin modelling of Lagha &; Manneville (2007) . 



A^y = 15 is found to be an acceptable compromise allowing the consideration of domains wide 
enough to accommodate several bands. 



2. Results 

We perform direct numerical simulations using Gibson's public domain code ChannelFlow 
(www.channelflow.org). This pseudo-spectral Fourier x Chebyshev x Fourier code is used at 



voluntarily degraded resolution, understood as a consistent modelling strategy (Manneville &; 
Rolland, 2011). Preliminary work showed that the band regime is remarkably robust, see Fig.[lj 
For systems with in-plane dimensions Lx (streamwise) and Lz (spanwise), the main features 
of transitional PCF are appropriately reproduced, including the pattern's wavelengths, with 
Nx — Lx and N^ — 3L^ collocation points, while taking Ny — 15 Chebyshev polynomials 
just yields a downward shift of [Reg,Ret] of about 15% to ^ [275,350]. Using this optimal 
compromise between accuracy and computational load, we consider a pretty wide domain of 
size Lx — 432/i and Lz — 256/i with the power of a desktop computer only. 

Starting from a state in the featureless regime at Re = 450 » Ret, the Reynolds number 
is decreased by steps down to the range of interest around Rcg '^ 275. At Re = 275 a steady 
pattern with 3 oblique bands is obtained with streamwise period A^; = 144/i and spanwise period 
\z — 85.3/i, in reasonable agreement with experimental findings by [Prigent et al. (2003). When 
Re is decreased further, the pattern breaks and decays. Figure [2] illustrates the different steps 
observed during the transient that starts from a state at Re = 275 upon setting Re = 272.5. 

Two different processes are seen to take place: the breaking of a continuous band and the slow 
retraction of band fragments at roughly constant speed. These processes are always present in 
the experiments that have been undertaken, some of which are featured by means of time series 
of the distance A to laminar flow in Figure |3j Band breaking is crucial in controlling the decay. 
For example, the transient at Re = 273.75 is marked by the persistence of the 3-band state for 
about 7000 /i/t/ before two bands break nearly simultaneously (green full line) and recede. In 
contrast, in the transient at Re = 272.50 pictured in Figure [2J bands break one after the other 




Figure 2. Snapshots taken during the decay at Re = 272.5. From left to right and top to 
bottom, t = 0: initial state, three continuous bands (1,2,3); labelling takes periodic boundary 
conditions into account, t = 2250: band 1 breaks, t = 5250: band 1 recedes, t = 6000: two 
continuous bands (2,3). t = 15750: band 2 breaks, t = 18750: band 2 recedes, t = 20250: one 
continuous band (3). t = 21750: band 3 breaks, t = 25500 & 26250, band 3 recedes and decays. 
Time is reset at the beginning of the experiment when Re is quenched from 275.0 to 272.5. The 
streamwise direction is along the vertical. Colour scale for local perturbation energy from blue 
for laminar to red for peaks of turbulence. Streamlines of the large scale fow, conspicuous in the 
laminar regions, are indicated. This flow is very weak except close to the turbulent patches. Its 
skewed quadrupolar shape around the tips of turbulent band fragments is typical. 



and the 2-band state remains unbroken during about 8000 h/U (red dashed line). The three 
bands break shortly after the beginning of the experiments at Re = 271.25 and 270 (brown and 
magenta lines). 

While band breaking leaves essentially no trace on the time series of A (Fig. [s]), the slow 
withdrawal of the band fragments corresponds to the linear decrease of this quantity which 
is roughly proportional to the cumulated length of the band fragments. The decrease is then 
easily understood as retraction at constant speed, itself a function of the number of fragments 
since retraction takes place at each end of a given fragment, statistically. This explains that the 
decay of turbulence is faster for Re = 273.75 with two simultaneously broken bands than for 
Re = 272.5 with only one band broken at a time. Both processes, the breaking of bands and 
the trimming of band fragments at their extremities are attributed to the chaotic dynamics at 
the MFU scale: Being transient in the range of Re of interest, it implies a finite probability for 
turbulence to decay locally. 

Band breaking happens when a turbulent band collapses in a domain of size at least as large as 
the band itself. This specific process has been studied in a smaller domain of size 144x84 able to 
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Figure 3. Time series of the distance to laminar flow during transients started from the same 
initial condition and different values of Re in the wide system (432 x 256). 
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Figure 4. Time series of the distance to laminar flow for the sustained banded state at Re = 275 
in the small system (144 x 84). 



contain just a single wavelength of the pattern (Rolland & Manneville, 2011). Figure [4] displays 
the time series of the distance to laminar flow for Re = 275 during an experiment of very long 
duration (3x10^ h/U) over which this specific turbulent regime remains sustained. Excursions 
toward low values of the distance, such as the one at t ^ 8500 h/U are not sufficiently dangerous 
nor sufficiently frequent to trigger the decay of the system to laminar flow. By contrast, for 
Re < 275, decay is observed in finite duration simulations, owing to the more frequent occurrence 
of more dangerous excursions. 

Probability density functions (PDF) of the instantaneous value of A are then constructed. 
They are displayed in Figure [5] for different values of Re, where it can be seen that the most 
probable value slowly decreases with Re but, more importantly, that the probability of low 
values, corresponding to the dangerous excursions, increases drastically in a narrow range of 
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Figure 5. PDFs of the distance to laminar flow in the small system for different values of Re. 



Re. (Irreversible decay takes place for distance < 0.13.) Decay is not observed for Re > 275 for 
time series longer than 3x 10^, though the exponential tail in the corresponding PDF suggest 
that this event could occur with very low probability, i.e. could be observable in much longer 
simulations only. 

After a sufficiently wide gap has been nucleated, the steady retraction of the turbulent 
domains is then easily understood as a noisy propagation problem for a stable state (laminar) 
that gains over a metastable one (turbulent). This process can be understood in terms of 
spatiotemporal intermittency below threshold ( jChate fc Manneville , 1987). Care should be taken 
that the band fragments are surrounded by a large scale mean flow, which somehow complicates 
the dynamics. 



3. Discussion 

A detailed understanding of turbulence decay in systems of practical interest asks for going 



beyond the temporal framework underlying the chaotic transients paradigm (Eckhardt et a/., 



2008). Although it is essential to notice that the bands are observed in an range of Reynolds 
numbers where chaos is transient, it seems necessary to include spatiotemporal effects. The 
breakdown of turbulent bands then appears to be governed by a stochastic process involving 
coupled units featuring small clusters of tightly coupled MFUs with chaotic dynamics. Chaos at 
such scales and its possible local collapse are at the origin of the large deviations able to open 
gaps in the turbulent bands and are also responsible for the subsequent slow withdrawal of band 
fragments. These observations substantiate Pomeau's conjecture (Pomeau, 1986) about the 



relevance of nucleation processes in first-order thermodynamic-like phase transitions on the one 



hand, and of spatiotemporal intermittency as a directed-percolation stochastic process (Kinzel, 



1983) on the other hand, while showing that the complete hydrodynamic problem keeps some 



specificities since the role of large scale flow corrections within laminar domains and the origin 
of the bands are so far not elucidated. 
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